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The self-diffusion constant D is expressed in terms of 
transitions among the local minima of the potential (inher- 
ent structures, IS) and their correlations. The formulae are 
evaluated and tested against simulation in the supercooled, 
unit-density Lennard- Jones liquid. The approximation of un- 
corrected IS-transition (1ST) vectors, Do, greatly exceeds D 
in the upper temperature range, but merges with simulation 
at reduced T ~ 0.50. Since uncorrelated 1ST are associated 
with a hopping mechanism, the condition D ~ Do provides 
a new way to identify the crossover to hopping. The results 
suggest that theories of diffusion in deeply supercooled liquids 
may be based on weakly correlated 1ST. 



I. INTRODUCTION 

The configuration r(t) of a liquid may Q be mapped to 
the local minimum R(i) , of the potential energy U (r) , to 
which it drains; both vectors denote the 3N coordinates 
of a system with N atoms. Stillinger named the minima 
inherent structures (IS) and Stillinger et. al. have ||] 
pioneered the application of this mapping to the theory 
of liquids. The configuration space decomposes into the 
basins of attraction of the IS, and so too does the con- 
figuration integral, providing a new approach to equilib- 
rium statistical mechanics. Turning to dynamics, Still- 
inger and Weber JlJ calculated the average T-dependent 
rate < 0Ji S (T) > of IS-transitions (1ST) for an atomic liq- 
uid, obtained an activation energy comparable to that for 
diffusion, and generally discussed the connection between 
1ST and physical dynamics. 

With the system described by both the continuously 
varying r(t) and the discontinuous R(i) an explicit real- 
ization of H Goldstein's ideas about dynamics over the 
potential energy landscape is obtained. The most chal- 
lenging area of liquid theory is the supercooled state, 
where the landscape paradigm is widely accepted. Below 
a critical temperature T c the system is believed to re- 
main in a basin for relatively long times, with infrequent 
activated hops to neighbor basins. At T > T c motion is 
freer and not activated. The character of dynamics has 
been related to {| the T-dependence of the IS energy, 
and to H the roughness of the landscape. Neverthe- 
less, no detailed calculation the self-diffusion constant D 
has been given. In this Letter we demonstrate the re- 
lation between D and 1ST in supercooled unit-density 
Lennard- Jones (LJ). At the lowest attainable T, D is de- 
scribed accurately by a simple Markov approximation, 
and a qualitative change in the dynamics is evident. 



II. DIFFUSION AND INHERENT STRUCTURE 
DYNAMICS 

The diffusion constant is determined by the long-time 
linear behavior of the mean-square displacement (MSD). 
Since the magnitude of the S return distance u(t) = 
r(t) — R(i) is bounded by the size of a single basin while 
the MSD increases without limit the IS-MSD also serves, 
< (AR(t)) 2 > /6N = Dt,t-> oo, where AR(t) = R(f)- 
R(0). The IS-MSD is the sum of the 1ST vectors, the 
separations of successive IS, AR(t) = X)a=i after 
n transitions. Squaring, averaging for fixed n, averaging 
over the distribution of n at time t and dividing by 6N 
yields 

oo 

D = [< {SK 2 )lu is >+2^< (<5R Q ■ SR a+f3 )uj ls >]/6N, 

13=1 

(1) 

where the sum is the correlation function C{0) for /3- 
neighbor 1ST vectors and we have assmed that t is much 
longer than the persistence time of any correlations, re- 
quired for diffusive behavior. Approximations to D keep- 
ing P < m only are denoted D m and the IS-Markov ap- 
proximation is Do, the first term on the RHS. The 1ST 
vector correlation C(/3) is somewhat analogous to the ve- 
locity correlation. 

With large times and numbers of trajectories the aver- 
ages in Eq [l] factorize into < uj is > times averages of <5R 
only, and in particular the Markov approximation yields 
D =< 5R 2 >< uj is > /6N. In computer simulation of 
supercooled liquids fluctuations are large and we use Eq|l| 
as is in averaging over multiple MD runs; nevertheless 
the difference is within the noise and our results may be 
understood via the factorization approximation. To see 
how D achieves its status as an intensive quantity con- 
sider that a liquid with short-ranged correlations may be 
H roughly divided into independent local regions. Suc- 
cessive IS differ Q] in the coordinates of a small number 
of particles and < <5R 2 >~ 0(1); 1ST are local rear- 
rangements. On the other hand the number of regions 
is O(N) and a change in any region changes the IS of 
the entire system, so |lj < u> is >^ O(N) and D ~ O(l). 
The N-dependence of < Ui S > indicates that IS-dynamics 
studies should use the smallest realistic system to avoid 
an intractable 1ST rate with no direct correspondence to 
a physical rate. 
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III. SIMULATION METHODS AND RESULTS 

We have previously studied jjj the unit-density super- 
cooled LJ liquid, and use it with N — 32, and with 
methods from prior work, to test the relation between 
IS-dynamics and diffusion. Natural LJ units are used 
throughout, well depth e for T, hard core radius a for 
distance, time unit tlj = (mcr 2 /e) 1 / 2 (2.18 ps for argon) 
and m is the mass; the crystal melts at T ~ 1.6. A hot 
liquid at T — 5.00 is cooled in one step to a temperature 
in the 1.20-0.60 range. The system is equilibrated for 2.5 
tlj, data are gathered for 62.5 tlj, T is decreased by 
0.02 and the process is repeated 10-25 times, generating 
a single quench run. The cooling rate is 3.8A10~ 4 . Dif- 
ferent T — 5.00 configurations lead to different behavior, 
and thus we obtain an ensemble of quenches. At N = 32 
the abrupt drops in U signaling solidification, common 
at N = 256, do not occur but some quenches develop 
solid-like pair distributions and these are discarded. Re- 
sults are averaged over all surviving quenches, typically 
5-15, at each T. Quench- to- quench fluctuations are much 
larger than any systematic changes over T = 0.02, so we 
also average results at each T with those from the next 
higher and lower T. Even so our data are somewhat noisy 
but the trends are clear for 1.10 > T > 0.34. For now 
T = 0.34 is the lower limit because of the scarcity of 
liquid-like states and the magnitude of fluctuations com- 
pared to the mean in dynamical quantities (D,u>i s ). 

Conjugate gradient (CG) minimizations are per- 
formed every 5 time steps (dt — .00125), or 160 
minimizations/TXj. Since the range of < uji S (T) > is 
from 9.6 IST/r LJ at T=1.10 to 0.28 at T=0.38, this 
should be sufficient. The determination of whether a 
transition has occurred is an important, nontrivial mat- 
ter. We begin, following Q Donati et. al., by calculating 
the distribution g(d) of distances d — log[((R(t) — R(t — 
5dt)) 2 /A) 1 / 2 ] between the current and prior IS, with no 
reference to the presence or absence of a transition. Dis- 
tributions for T = 1.10, 0.90, 0.60, and 0.40 are shown 
in Figure |. There is a large peak, not shown, 
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FIG. 1. Distributions of logarithms of IS-transition dis- 
tances. Top to bottom at d = -1,T = 0.90,0.60,0.50,0.40. 

around d = —6 which identifies the numerical uncer- 
tainty of the algorithms when the IS has not changed. 
Transition distances exhibit a bimodal structure, with a 
small peak at d ~ —4.5. The spikes in the larger, low-d 
lobe are not noise but represent specific, frequently occur- 
ring separations. They vanish abruptly between T = 0.60 
and T = 0.40, the range in which we believe the crossover 
to hopping begins. 

Not all 1ST are associated with diffusion. A change in 
state of a two-level system will change the IS. Stillinger 
has argued Q that motion among megabasins is required 
for diffusion, with motion among the basins comprising 
a megabasin non-diffusive. Transition vectors with large 
d are most likely to reflect diffusion. Similar considera- 
tions have H entered the efforts to relate D to Im — uj 
instantaneous normal modes, where non-diffusive modes 
must be discarded. Accordingly, we count an 1ST when d 
falls in the high-ei lobe of g(d), specifically for d > —1.55, 
marked with an arrow in Figure [l| This choice has little 
impact on our estimate of D, because increasing < LOi S > 
by taking more small-d 1ST would decrease < (^R 2 ) >, 
leaving the product almost unchanged. 

At each T we obtain several 1ST quantities, and D de- 
termined as usual from the MSD. Figure |^ compares sim- 
ulated D with the Markov approximation Do and with 
Di, including correlation of adjacent 1ST vectors. At 
T ~ 1 Do >> D, direct evidence of 
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FIG. 2. Diffusion constant vs. T. Top to bottom at T = 1, 
Do (IS-Markov approximation), D\ (adjacent 1ST correla- 
tions), simulation. 

large negative correlation. However adjacent correla- 
tions only produce a ~ 10% correction, so C{(3) must be 
'long ranged'. 

As T decreases the situation changes strikingly. While 
Dq(T) initially decreases much faster than D(T) and 
would extrapolate to zero at T ~ 0.6, it abruptly changes 
slope starting at T ~ 0.7 and begins to merge with 
D at T ~ 0.5 (we believe the bump at T ~ 0.6 will 
vanish with more averaging). This crossover reflects 
behavior of < w, s >; < (<5R 2 ) > is roughly linear 
for 1.10 > T > 0.34. Figure | highlights the region 
0.54 >T> 0.34. The joining of the curves at T ~ 0.5 is 
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FIG. 3. Low-T region of Figure 
tion improves with decreasing T. 



IS-Markov approxima- 



clear, and D (T) and D\(T) are approaching quan- 
titative accuracy below T ~ 0.4. In the same T-range 
< Ui S (T) >, Figure [|, is undergoing its abrupt fall to the 
bottom of the 
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FIG. 4. Averaged IS energy vs T, data and Bezier fit. 

landscape, previously B associated with crossover to 
hopping. Figure |4| is similar to that given JhJ by An- 
gelani et. al. for unit-density 'modified LJ', N = 256, 
where they estimate T c = 0.475. 



IV. DISCUSSION 

In a hopping mechanism the system is constrained by 
the need for activated barrier crossing to spend a long 
time (e.g. compared to a vibrational period) in a basin 
between hops; it is plausible that they should lose cor- 
relation during the wait with Dq ~ D. On the other 
hand at higher T, with thermal energy comparable to 
barrier heights, 1ST are 'bookkeeping' events as the sys- 
tem moves freely across IS-boundaries. As Stillinger and 
Weber jlj discussed, a burst of 1ST can then be gener- 
ated by a small displacement through a region of closely 
spaced boundaries. If the 1ST rate is unpyhsically large 
accordingly, so too will be D , and negative 1ST- vector 
correlations C{j3) will be required to obtain the correct 
D. 

Our results are consistent with the above scenario and 
suggest that the transition to a hopping mechanism oc- 
curs over the range 0.50 > T > 0.34, possibly associ- 
ated with the fall of the system to lower- lying IS. The 
crossover is usually identified by fitting D(T) to the pd| ] 
mode-coupling form, but the criterion D ~ Do is perhaps 
a more direct indicator. We cannot see characteristic Q 
deeply-supercooled behavior in this simulation of a sim- 
ple liquid. Irregardless of any interpretation, it is very 
encouraging that Dq becomes a good first approximation 
at the lowest available T. Of course correlations also ex- 
ist at low T. If |||l2) the syste m revisits a small group 
of basins over and over, perhaps due to low connectivity, 
the corresponding 1ST will not cause diffusion and C (/3) 
must be negative. The hope is that, for a broad spectrum 
of materials, such physical effects will be more tractable 
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than the gross overestimate of D at high T. 
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